The Lund Model at Nonzero Impact 

Parameter* 



CO 

o 

CN ; Romuald A. Janik a , Robi Peschanski 6 ^ 
-t— > 
o 

a Jagellonian University, 

00 ' 

' Reymonta 4, 30-059 Krakow, Poland. 

^ i b CEA/DSM/SPhT Saclay 

^* . Unite de Recherche associee au CNRS 

O ■ CEA-Saclay, F-91191 Gif/Yvette Cedex, France. 

O : February 1, 2008 

co 
o 

a,: 



Abstract 



We extend the formulation of the longitudinal 1+1 dimensional 
Lund model to nonzero impact parameter using the minimal area as- 
sumption. Complete formulae for the string breaking probability and 
^ ' the momenta of the produced mesons are derived using the string 

■ worldsheet Minkowskian helicoid geometry. For strings stretched into 

the transverse dimension, we find probability distribution with slope 
linear in similar to the statistical models but without any ther- 
malization assumptions. 



1 Introduction 

One phenomenological ingredient of models for hadronic reactions is the de- 
scription of the hadronization process when the (hard) partons produced at 
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the interaction point get dressed to form the observed hadronic particles 
which are registered in the detector. This is a theoretically difficult prob- 
lem due to its nonperturbative character, however phenomenological models 
are quite successful in the description of data. Models which describe this 
behaviour are essential for the comparision of QCD predictions for partonic 
processes with the experimental data. 

One of the most successful hadronization models is the Lund model ^ 
12] based on the effective string picture of QCD at strong coupling. The 
conventional Lund model 1 gives a description of the break-up of a colour 
string into produced hadronic states. It is directly applicable in the context 
of processes where the initial qq pair is produced at a single point, the partons 
then fly away from each other, a colour string is stretched between them, 
which eventually breaks up producing hadrons. A fundamental property of 
the confining strings is the so-called area law jl] , from which the probability 
formula in the Lund model may be advocated. 

The model is basically 1+1 dimensional, transverse momenta of the gen- 
erated mesons are taken into account by substituting the masses m 2 by 
m\, = m 2 + k\. In particular the transverse momenta are considered to be 
in principle randomly distributed, decorrelated along the chain of mesonic 
emissions and also uncorrelated with rapidity. 

The aim of this paper is to show that the whole 1+1 dimensional setup 
of the Lund fragmentation model based on the area law can be generalized 
to the 3+1 case when the initial quark and antiquark are also separated by 
a transverse distance L (impact parameter). These types of colour string 
configurations may occur between coloured remnants of a projectile and of 
the target. They could also appear locally in the hadronization of a partonic 
cascade between neighbouring partons at some transverse distance. 

Our results are the following. We give the full mathematical derivation of 
the probability distribution function for the production of mesons with their 
4- momenta (([11)1 and formulae (|22 p -([24| ) ). This is based on the Minkowskian 
helicoid geometry which is the relevant string worldsheet for nonzero impact 
parameter L. The results depend on the key parameter p = x/L, where x 
is the rapidity interval. For large p we recover the Lund model probability 
distribution, while for small p new features appear. We obtain a linear slope 
in irtT and not m\ as is usually assumed in effective string models. This 

1 By the 'conventional' Lund model we mean the original 1+1 dimensional formulation 
not containing the extension to gluonic strings 3 . 
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derivation based just on classical considerations has to be contrasted with 
similar predictions based on statistical model assumptions. 

The plan of this paper is as follows. We will first review the standard 1+1 
Lund model along the lines of [3]. Then in section 3 we will introduce the 
Minkowskian helicoid minimal surface which is crucial for the construction 
of the generalization to nonzero impact parameter. In section 4 we will show 
how the basic formulae of the Lund model get modified in the new setting. 
In section 5 we will analyze the limiting cases of small and large p. In section 
6 we close the paper with considerations of phenomenological interest. 



2 1+1 dimensional Lund model 

Suppose that a highly energetic qq pair is produced at the point x^ = 0. 
These particles move along straight lines 

x = yto x = -y^o (1) 

with relative rapidity 

Y = — log . 2 

A 2 l-v K ' 

According to the effective string picture of QCD, a colour string is formed 
between them. Since by assumption we are in the confining regime, the 
string lies in the most energetically favourable configuration i.e. it forms the 
minimal surface joining the two lines (see fig. la). This surface is of course 
flat in this case. This follows from the Nambu-Goto action for a confining 
string worldsheet Sng = K • Area. 

The process of hadronization is described by the break up of this string 
(see fig. lb) forming a pair of massless quarks, each of which moves along 
the light cone directions and which get converted into hadrons (at the upper 
tips in fig. lb). The probability of such a breakup configuration is given by 
the area law: 

Probability ~ exp (—6 ■ Area) (3) 

where Area is the area spanned by the string and b is a phenomenological 
parameter which is related to the break-up probability of the string per unit 
area 0. 

The area can be calculated as a sum of the areas of the diamond-shaped 
sectors (like the one in fig. lb between the dotted lines). Denoting the 
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Figure 1: Configuration with zero impact parameter, a) Initial qq configura- 
tion; b) a sample configuration of the confining string with break-up points 
(the dotted lines divide the string area into diamond-like sectors - see text). 



coordinates of the adjacent points of breakup by x*) and (tf,Xf) one 
obtains the formula 

Area if = J (2x^x + f -x +f x- f - x +i x_ f ) (4) 

where 

x± — t ± x . (5) 

At this stage the essential probabilistic set-up of the Lund model is com- 
plete. 

The second ingredient in the Lund model is a prescription how to as- 
sociate the produced physical mesons to the string break-up configuration 
characterized by the break-up points (tk,Xk). 

Let the string worldsheet be parameterized by the coordinates r and a. 
Then the Nambu-Goto lagrangian 



C = K^j- (h TT h aa - h? Ta ) (6) 

is expressed in terms of the induced metric 

Kb = daX^dbX^ . (7) 

The 4-momentum of the produced meson formed between the adjacent 
points of breakup (tj,Xj) and (tf,Xf) follows from the formula for the mo- 
mentum of a (classical) string stretched between these points [6J: 

P m = / 2±L da _ ^L dT ( 8 ) 
Jcdx^ dx"" K ' 
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where C is (any) curve joining the breakup points. For flat {h ah = r] ab ) or 
more generally conformally flat coordinates (h a b = e^^'^rjab) this simplifies 

to 

P» = K / (x^da + x'^dr) . (9) 
Jc 

Applying (jHJ) to the flat minimal surface shown in fig lb yields 



,7 ; 



P U = K( Xf 
P 1 = K(tf-U) 

P ± = 0. (10) 

The final ingredient, which is not derived from the model, is the mass shell 
condition for the produced mesons. Thus 5(P 2 — m 2 ) is put in where m 2 are 
the physical meson masses. 

The resulting probability distribution function PQ for the production of 
mesons with momenta {pj} is thus: 

dP n ({ Pj }, P tot ) ex d 2 Pj 5(p 2 - m 2 ) 6 \^2 Pj - P t0 )j exp(-b ■ Area) (11) 

where 

Area = ^ Areaij (12) 

is the total area of the string worldsheet (see fig. lb). 

Note that all the above steps and formulae, apart from the mass-shell 
condition, follow just from the geometry of the flat surface between the ini- 
tial lines, which in turn is unambigously determined as the minimal surface 
between the trajectories of the initial energetic qq pair which generates the 
initial colour string. The point that we would like to emphasize here is that 
it is possible to use this property to generalize the model to configurations 
with nonzero impact parameter. 

3 The minimal surface geometry: the Mink- 
owskian helicoid 

When the string is formed between a quark and antiquark, moving with 
relative rapidity \, and separated by a transverse distance L, the string 
worldsheet will stretch along the minimal surface formed by the two lines. 
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Figure 2: A string breaking on the helicoid. a) picture projected onto the 
longitudinal plane (to be compared with fig. lb); b) the same configuration 
in the r — a variables. 

An Euclidean version of this surface has been used in an approach to study 
high energy scattering using the AdS / CFT correspondence [7] . 

The minimal surface here 2 is a (Minkowskian) M-helicoid which can be 
parametrized by 

t = - sinhpr cosh pa 
V 

x = - sinh pr sinh pa 
p 

x± = a (13) 

where p = x/ L, a = — L/2 . . . L/2 and r starting 3 at 0. The string at t = 0, 
when the particles are separated by a transverse distance equal to the impact 
parameter, is stretched along a straight line. 

The geometry of the helicoid is encoded in the induced conformally flat 
metric: 

h ah = (cosh 2 pr) • r, ab . (14) 

2 The assumption that the string at t = is a straight segment in impact parameter 
selects the helicoid minimal surface. If the initial string would be in some excited con- 
figuration the analysis would possibly lead to different solutions. Note that an a priori 
different transverse configuration, going beyond the original Lund model, has been studied 
in [H] with a numerical treatment based on the ARIADNE 2 program j^]. 

3 Our formalism can be easily extended to other initial conditions for r. 
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Hence 



ds 2 = cosh 2 pr • (dr 2 — da 2 ) 



(15) 



Since the produced quarks are assumed to be massless, they follow light-like 
trajectories (ds 2 = 0) on the helicoid which are just straight lines at an angle 
of ±45° in the r — a variables. The area element is 

dA = V— det h = cosh 2 pr drda . (16) 

It is also convenient to look at the helicoid projected onto the longitudinal 
t — x plane. Then one can show that 

— = tanhper pa = - log — . (17) 

t 2 x_ 

So the lines of constant a correspond to fixed longitudinal velocities 4 v. 

The t coordinate can also be related to the longitudinal coordinates 
through 

sinhpr = p^/x^xZ . (18) 

As we saw in the previous section all ingredients of the Lund model can be 
formulated as purely geometric constructions based on the minimal surface 
which represents the colour string. We thus have to redo the calculations of 
section 2 substituting the helicoid for the original flat surface. 



4 The Lund model at nonzero impact param- 
eter 

We will now derive the formulae for the break-up probability and produced 
momenta for the nonzero impact parameter case which are enough to define 
the generalization of the full probability distribution function pi)). The only 
thing that changes is the underlying geometry which is no longer flat but 
is given by the Minkowskian helicoid. It is convenient to label the breakup 
points by the longitudinal coordinates (t, x) (i.e. to look at the helicoid 
projected onto the t — x plane) and subdivide the area into same sectors as 
before (the boundaries being lines of constant slope x/t = v = const, i.e. 

4 This is also the definition of the boundaries of the sectors in fig. lb in the original 
zero impact parameter model. 
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constant a) . In fig. 2a we show the breakup in the t — x plane while in fig. 2b 
we depict the same process 5 . . 

We calculate first the area of the sector defined by the breaking up points 
(er,j, Ti) and (07, r/). The coordinates of the tip of the sector are: 

^ = n = Tf + n + °f-°i (19) 

The area is then given by 

f a * fTi+a-Oi rcr f rT f -cr+cr f 

Areaif = da cosh pr dr + da cosh pr dr (20) 

J(Ti JO Jcr, JO 

Using the integral 

f s , 9 , s sinh(2ps) ,^ , 

/ cosh 2 prrfr = - + ^—i- (21) 

the Areaif can be calculated in closed form: 

Area^ = - ((a/ - a { ) 2 - (r/ - t^) 2 + 2(a/ - <7 f )(r/ + n)) + 

+ — (2 cosh(p(a/ - <Ji + Tf + n)) - cosh(2pr/) - cosh(2prj)) (22) 

We now have to compute the momentum carried by the string between 
the breakup points. This is identified with the 4-momentum of the produced 
meson. Using formula we get 

P° = — (cosh pTf sinh pa / — coshpTj smh pai) 
P 1 = — (cosh pTf cosh pa f — coshprj coshpcrj) 

P± = K{r f - n ) . (23) 



5 We should emphasize that in contast to the break-up points, which have a definite 
physical meaning, the subdivision of the string worldsheet into sectors is to a large extent 
arbitrary. Our choice here was to subdivide the area in a way similar to the flat string case 
in order to facilitate the comparison with the 1+1 dimensional Lund model. Of course the 
final probabilistic formulas contain only the sum of all sectors which is uniquely defined 
in the r — a variables. 
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It is convenient, also for comparison with the zero impact parameter case to 
rewrite this in terms of longitudinal coordinates 



P 1 



K 

K 

K 
P 



1 + 



1 



1 + 



1 



p^x +i x. 



^ ~ ' p 2 x +f x^ f f V 

(log [y/l +p 2 X +f X- f + Py/X +f X- f 



-U 



(/ 



(24) 



Formulae (}2"2*|) . (|2*Hj) and equivalently (|24|) when inserted into (fTTj) com- 
pletely define our model and can be used as a basis of a Monte-Carlo simu- 
lation. 



5 The limits of small and large aspect ratio p 

The key parameter which characterizes the surface geometry is p = \l L. The 
scale of the momenta of the produced mesons is set by the string tension k. 

Small impact parameter 

Let us first verify that when p — > oo (large rapidity or zero impact parameter 
limit) we recover the formulae of the conventional Lund model. Indeed it 
is clear from (|24j) that the momenta of the particle produced go over to the 
conventional Lund expressions (note that P± ~ (logp)/p — > 0). In fact we 
have to impose the condition that x± of the break-up points are kept fixed 
in the limit. From (j!7j) this means that o —>■ with pa fixed. 

The Area of each sector (|2*2*j) also reduces to the Lund expression (|3J). To 
see this we note that for large p (and keeping the longitudinal coordinates 
of the break up points fixed) the polynomial terms in (|2*2*|) vanish, the only 
contribution comes form the hyperbolic cosines. This can be easily evaluated 
using e pa = \Jx + /x~ and e pr ~ 2p^/x + X- (for p — > oo). The result is the 
Lund expression (jlj). 

Since we obtained the expression (jlj) where there is no explicit remaining 
p dependence, our assumption that the break up longitudinal coordinates 
do not scale with p but remain of order 1 is self-consistent. In this way we 
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verified that the conventional Lund model is a smooth limit of the nonzero 
impact parameter generalization. 



Large impact parameter 



It is especially interesting to study which new phenomenae may arise in 
comparision to the original model, and whether these are phenomenologically 
important. To answer this question in full one has to perform a Monte Carlo 
simulation. Here we will just consider the limit of large impact parameter 
(small p limit) where one can expect the largest deviations from the original 
Lund model. 

Before we make the analysis let us note that the string worldsheet in the 
p — ■> limit is a flat string stretched in the transverse plane which is still 
unbroken at t = 0. Then it will fragment into mesons which will have nonzero 
transverse momenta even at the classical level. 

Let us check how does the area scale when p — > 0. If we where to keep the 
longitudinal coordinates fixed and of order 1 we would arrive at a singular 
area 6 ('action') for p — > 0. Hence the above assumption in this case is 
unjustified. The correct choice is in fact that the coordinates (a, r) remain 
of order 1. This is natural physically as the transverse mass in the small p 
limit stays finite and becomes my = kA<j = — Cj). The r coordinate is 
also directly linked to the transverse momentum through At = P±/ K. 

Considering now a and r fixed and starting with the area formula 
one gets for p — > 0: 



Area 



if 



1 r 



AV - A 2 r 



t ' Act ~ - — - + r 



m 
4^ 2 



1 

K 



(25) 



where 



Ac 



A 7 



(26) 



Therefore the break-up probability for the formation of a meson takes the 



3 One gets then 



Area x , 



x -i x +f + x - f x +i 



i/ x -i x +iy /x - f x +f 



32p 2 



log 



x +f 



log 
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following simple form 



Probability oc exp ^— ' ex P ^ ~~ T+ ) ' 

where r + can be related to the mean time value of two adjacent break-up 
points in the limit p — > 0. Indeed, one has 



tf — xf + yt 2 f — xj- = - sinhprj H — sinhpr/ — > Ti + r/ w/ien p — > . (28) 

6 Phenomenological considerations 

Formula (j2*Tj) has quite intriguing properties. Let us assume that the value 
of r + fluctuates around some average value (r). We emphasize that this 
assumption should be checked in a complete Monte-Carlo simulation. Then 
the distribution of produced particles is predicted to follow a simple formula 

with a slope linear in m^. This is quite unusual in string based models where 
generically a behaviour quadratic in m% appears 7 JU] i n relation with the 
Schwinger mechanism [llj. Here the slope linear ttit arises due to subtle 
cancellations between quadratic terms in the helicoid area and thus appears 
directly at the classical level. 

It is interesting to compare these properties with those of the generic sta- 
tistical model ^2] which successfully describes both the slope in rriT and the 
abundances of different particle species in high energy reactions |13j . In this 
model the creation of particles is supposed to come from thermalized clus- 
ters of limited center of mass energy. Our model, at small p, leads to similar 
conclusions but without invoking the thermalization assumption. Indeed we 
have checked that the slope in itlt and the average relative distribution of 
the multiplicities of different species can be qualitatively described by a sim- 
ple formula of the form (|29jl. In order to go further and to determine the 
parameters like (r) one needs to perform a complete Monte-Carlo simulation 
which goes beyond the scope of this paper. 

7 In ref. jlOj it was suggested that the phenomenologically preferred behaviour linear 
in m-T can be obtained by considering gaussian fluctuations of the string tension. 
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Let us finally note another possible implication of the string breaking 
model based on the helicoid geometry. Since the helicoid is extended in the 
transverse direction, transverse momenta get produced at the classical level 
on the same footing as longitudinal ones. This opens up an exciting possibil- 
ity that correlations between p t and rapidity could be naturally incorporated 
in the model. We leave a detailed investigation of these issues which should 
be based on a Monte Carlo simulation for future work 8 . 
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